Detecting the tunneling rates for strongly interacting fermions on optical lattices 
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Strongly interacting fermionic atoms on optical lattices are studied through a Hubbard-like model 
Hamiltonian, in which tunneling rates of atoms and molecules between neighboring sites are assumed 
to be different. In the limit of large on-site repulsion U, the model is shown to reproduce the t-J 
Hamiltonian, in which the J coefficient of the Heisenberg term depends on the particle-assisted 
tunneling rate g: explicitly, J = 4g 2 /U. At half-filling, g drives a crossover from a Brinkman-Rice 
paramagnetic insulator of fully localized atoms (g — 0) to the antiferromagnetic Mott insulator of 
the standard Hubbard case (g = t). This is observed already in the number of doubly occupied sites 
under the intermediate coupling regime, thus providing a criterion for extracting from measurements 
the effective value of g. 
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I. INTRODUCTION 

In recent years, the field of ultracold atoms has proven 
to be a very rich and growing research area [1], one of the 
basic motivations residing in the possibility of experimen- 
tal access to and probing of a wide range of prominent 
problems in condensed- matter physics. For example, the 
rapid development of experiments on ultracold atomic 
gases loaded onto optical lattices has supported the ob- 
servation of the superfluid-Mott-insulator transition in 
the Bose-Hubbard model [2|. 

One of the most relevant challenges in this field turns 
out to be the simulation of the fermionic Hubbard model 
0, 13 an d its extensions Q , which are believed to be 
the major candidates in describing high-temperature su- 
perconductivity. Most recently, in very intriguing exper- 
iments with fermionic 6 Li [7] and 40 K [8| atoms loaded 
onto optical lattices, a Mott-insulating phase was found 
under the regime of strong on-site repulsion between par- 
ticles, characterized by the presence of one atom per site 
and by a number of doubly occupied sites (doublons), de- 
pending on the strength of the external parameters. Such 
a state turns out to be relatively easy to achieve experi- 
mentally -due to its incompressibility- and its observed 
behavior is qualitatively consistent with well-established 
properties of the Hubbard Hamiltonian [9j. In addition, 
it has been pointed out that under a regime of strong in- 
teraction, the necessary Feshbach resonance could induce 
highly nontrivial processes, which involve multiband pop- 



ulations and off-site interactions. These are described by 
a generalized Hubbard Hamiltonian proposed by Duan 
[10| , previously known in the literature as the Simon- 
Aligia model [11]. In fact, this model differs form the 
standard Hubbard model in the effective tunneling rates 
between neighboring sites, which are assumed to be de- 
pendent on the type of particles involved, at variance 
with the standard case. The fact that in experiments it is 
possible to control separately the dynamics of molecules 
(doublons) and that of single atoms is explicitly included 
in the Hamiltonian. Here we show that this fact can be 
used to induce a strong coupling regime already at mod- 
erate repulsive interaction U . Such a regime can possibly 
be observed in experiments, for instance, as a lower num- 
ber of doublons in the insulating phase, with respect to 
the standard Hubbard case. 

In the following, we first investigate analytically the 
strong-coupling regime for the Duan Hamiltonian at ar- 
bitrary tunneling rates and fillings. Then we focus on the 
half- filled case and explore -also numerically by means of 
DMRG simulations- the behavior of the insulating state. 
The ultimate aim is to infer from experimental data the 
effective value of the tunneling rates in real systems. 



II. LATTICE MODEL HAMILTONIAN 



The model Hamiltonian reads 



H = - ^ [t + Sg(n id . + rija) + Stn^n^} c\ a 
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where a = ±1 (a = —a) identifies the two internal states on a d- dimensional regular lattice and rii = + 
of the fermionic atoms, (ij) denotes two neighboring sites 
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U > U c (g) > 0. 

In the following, we investigate how the actual Sg value 
also affects the behavior of the system in the intermediate 
to strong coupling region U > U c . In fact, it will turn out 
that the choice of St is irrelevant to the strong coupling 
behavior. 



III. THE STRONG-COUPLING LIMIT 
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Figure 1: Processes between neighboring sites described by 
the three independent tunneling rates t, g, and t a d- Here 
slashed and filled circles represent single atoms (of arbitrary 
spin orientation) and for diatomic molecules (doublons), re- 
spectively. 



Here t describes the direct hopping of atoms of a given 
population between neighboring sites, while Sg = g — t 
and St = t + t a d — %g are the deviation from the direct 
hopping case (in which Sg = = Si), induced by corre- 
lations in proximity of a wide Feshbach resonance. More 
precisely, as shown schematically in Fig. [TJ g describes 
the tunneling configuration in which one atom is trans- 
ferred to a site already occupied by an atom of a different 
species; t a a accounts for the motion of one atom between 
two already-occupied sites, thus exchanging a molecule 
and a fermionic atom located at neighboring sites. Fi- 
nally, U is the energy cost of the molecule, which works 
as an effective detuning parameter. 

The derivation of Hamiltonian (pQ) for a system of ul- 
tracold fermionic atoms trapped on an optical lattice has 
been given in [10] and [12|. In the case of a double- well 
lattice, the numerical solution for two fermions interact- 
ing across a Feshbach resonance was also used to infer 
the range of parameters in which H works as an effective 
Hamiltonian [13|. 

In one dimension, the ground state of H has been in- 
vestigated in detail in the context of high T c supercon- 
ducting materials for the choice St = (see, e.g., Ref. 
P3] and references therein). Some results have also been 
achieved at St ^ pjHIll- Most recentl y t he case of 
imbalanced atoms has also been explored [T(J G3 . The 
ground-state phase diagram in the moderate interaction 
regime U < 4t differs substantially from what one would 
get by considering the standard Hubbard model. In par- 
ticular, for Sg > t/2 and U < U c (g), it is characterized 
by a phase that is superconducting even at half-filling 
and that displays an enhanced number of doublons per 
site, rid- By increasing the population imbalance, rid re- 
mains unchanged up to a critical polarization p c , beyond 
which some pair breaks and macroscopic phase segrega- 
tion takes place [l8|. At half- filling, the presence of such 
a superconducting phase shifts the insulating behavior 
characteristic of the Hubbard model to higher U values, 



To perform the strong-coupling limit of Eq. (pQ), we 
generalize to the present case a procedure developed in 
Ref. [19] for the Hubbard model. As the first step, we 
rewrite Hamiltonian (pQ) in terms of Hubbard projectors 
= where |a)j are the four states spanning 

the basis of the vector space at a given site j, with a — 
{0, +, -, 2} and |2)j = | + Explicitly, 
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We observe that the hopping terms in Eq. (|2]) can be clas- 
sified according to the change in the number of fermion 
pairs involved: 
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More precisely, To leaves unchanged the number of local 
pairs whereas T +1 (_ X ) creates (destroys) a pair. It is im- 
portant to note that = T_ m and [A^,T m ] = rriT rn , 
where Nd = ^Z^X? 2 , its expectation value being the 
total number of doublons in the system. In terms of 
the above operators, the Hamiltonian ((2]) reads: H = 
UNa + To+gp^+T-!). 

In the strongly interacting regime t/U << 1 we may 
decompose the fermionic Hilbert space of the model into 
subspaces characterized by different values of Nd'. 



N/2 
N d =0 



(4) 



Indeed, in this case the energy spectrum of H splits 
into well-separated subbands with different Nd, and a 
gap (~ U) opens between them. Hence, within this limit, 
we consider an effective Hamiltonian which does not mix 
different sectors of the Hilbert space, namely, conserving 
the total number of fermionic pairs. To obtain such a 
Hamiltonian we perform a unitary transformation that 
eliminates to the lowest order in U (the zeroth order) the 
terms not commuting with Nd. Explicitly, 



5 = 5( 1 > = ^(T+i-T_i), 



U 



(5) 
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so that, to the first order in U 1 , the rotated Hamiltonian 
reads 

= e s He~ s 

[S,H] [S,[S,H]] 

= ^ + ^r + 2! + • 

= /77V d + T + /7- 1 (^ 2 [T +1 ,T_ 1 ] 

+<7[r ,r_i] + <7[r + i,T ]) . 

Now one can think to iterate this procedure, eliminating 
the terms not conserving Nd also to the first order in 
Z7 -1 . This is implemented by choosing 

=gU- 2 ([T +1 ,T }-h.c.) , (6) 

so that the rotated Hamiltonian, up to second order in 
Z7 -1 , reads 

if ( 2 ) =UN d + T -^-[T- 1 ,T +1 ] + O(U- 2 ) . (7) 

The effective Hamiltonian in the lowest-energy sector 
(i.e., Nd = 0) is hence obtained as 

H eff = T (0) -JJ2 (SiSj - \n inj \ + 0(U~ 2 ) , (8) 

<ij> ^ ' 

where -as usual- the three-site term has been 
neglected [20]. Remarkably, H e ff is still in the form of the 
t-J model, as in the strong-coupling limit of the standard 
Hubbard case (g = t). However, in the present case the J 
coefficient has changed to J = 4g 2 /U. Note that by spe- 
cializing to half-filling one recovers the antiferromagnetic 
Heisenberg model. 

IV. DISCUSSION 

Equation dSj) provides -in the strong-coupling limit- 
the ground-state energy of the model Hamiltonian (pp) 
in arbitrary dimension, at arbitrary g and t a d starting 
from the known results for the t-J (generic filling) and 
Heisenberg (half-filling) models. Given the dependence 
of J on 5, Eq. (|5J) in fact describes the crossover from 
an antiferromagnetic insulator to a Brinkman-Rice [2l| 
paramagnetic insulator of fully localized fermions (at g = 
0). The effect was already predicted [17j in one dimension 
by means of numerical investigation of the model of Eq. 
(PQ), and here it finds its analytic confirmation. 

Such a result can be viewed as a consequence of the 
interplay of the exchange term (with coefficient oc g 2 /U) 
and the kinetic term describing the motion of single 
atoms in a background of empty sites. The latter, which 
is also known as the infinite- U Hubbard model, in one 
dimension describes the physics of a system of spinless 
fermions, with a vanishing number of doubly occupied 
sites, and energy per site e = — 2/7rsin7rn. In particular, 
at n = 1 it describes an insulator in which each fermion 
is fully localized at a different lattice site -no matter its 
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Figure 2: (Color online) Ground-state energy density e gs at 
half- filling, t = 1, versus U at three different g values, as 
obtained by numerical simulations with DMRG on L — 80 -r- 
120 sites. Dotted, dashed, and dot-dashed lines represent the 
results obtained for Hamiltonian (1), whereas solid lines are 
the results obtained for the Heisenberg model at the same 
parameter values. Asterisks denote in each case the value of 
U (U*) at which the two energies differ by less than 1%. Inset: 
e gs versus J at t a d — 0.8 (blue diamonds), t a d = 0.4 (green 
squares), and t a d = —0.8 (red triangles). 



spin orientation: the Brinkman-Rice insulator. Increas- 
ing g (and hence J) amounts to a gain energy, thanks to 
magnetic interactions, without affecting e: the insulator 
progressively increases its magnetic ordering and dimin- 
ishes its energy, to reach the Hubbard limit value. Even 
for Nd 7^ 0, inspection of Eq. still shows that higher g 
values correspond to lower energies, so that one expects 
the spinless character of the insulator to emerge (and Nd 
to diminish) with decreasing g. 

Some implications of the strong-coupling limit, ([?)) and 
dHJ), on the behavior of the system described by Eq. (pQ) 
at intermediate values of the repulsive interaction U are 
given here, with the help of numerical analysis. Since J 
also depends on g, one may expect that, with decreasing 
g, the strong-coupling regime -identified as the regime 
characterized by Hamiltonian (|5])- is reached at lower 
values of U. This is shown in Fig. [2j where we report 
the ground-state energy for Hamiltonian (pQ), compared 
to that of the corresponding t-J limit, Eq. flS}, at differ- 
ent g values and t = 1 . If we define [/* as the value of the 
interaction at which the difference between the two ener- 
gies is lower than 1%, it is seen that in fact U* decreases 
with decreasing g. This is in agreement with exact results 
on the two limiting cases known in the literature: while 
in the g = limit [16] the Brinkman-Rice insulator is the 
exact ground state at large enough finite U > U c (in one 
dimension, U c = 4t), for g = t the Heisenberg antiferro- 
magnet limit is reached exactly only at U = oc. From 
an experimental point of view, this observation suggests 
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Figure 3: (Color online) DMRG results for the number of 
doubly occupied sites (doublons) as a function of g at different 
tad values, with U = At. The solid horizontal lines represent 
the range of values of g consistent with a given value of rid- 
Dashed curves are guides for the eyes. 



ate coupling region U < U c (g) [17] where, as mentioned, 
rid could even increase slightly with decreasing g. 

Since rid is typically measured in experiments, Fig. [3] 
offers a straightforward way to extract information about 
the actual value of g in the experimental setup. We high- 
light with horizontal straight lines the region in which the 
range of values of g is compatible with a measured value 
of rid at a given U. Indeed, the range \t a d\ > t can be 
mapped into that with \t a d\ < t shown in figure. It is 
seen that at fixed U the standard Hubbard case reaches 
the highest relative values of rid- 

Besides, lower values of rid (and hence of g) are con- 
sistent with experimental data in [8] at moderate values 
of the energy trap, assuming that the dimension of the 
lattice enters only through the bandwidth. Note that 
already in the case of a double well [13], the value of g 
consistent with the observed spectrum was found to be 
lower than 1 (g « O.St). Moreover, in the case of real 
condensed matter systems, g « OM [9[. 



that the strong coupling regime could be explored even 
at a not too large detuning U, in the case where g is low 
enough. 

Moreover, since J in Eq. (jSJ) does not depend on t a d the 
same should happen for the ground-state energy versus 
J, at large enough U values. This is shown in the inset in 
Fig. [2j where in one dimension the ground-state energy 
is plotted and compared to that of the Heisenberg model 
at different t a d values and fixed U. It is seen that already 
at a moderate value of the repulsive interaction U, the 
effect of t a d is in fact negligible. 

In Fig. [3] the role of g for intermediate values of the in- 
teraction U is exploited to investigate the dependence of 
the number of doublons rid = Na/L on U. This analysis 
confirms that the strong coupling behavior also charac- 
terizes this regime, since rid still decreases by decreasing 
g, approaching smoothly the strong coupling value for 
g = 0. This result is consistent with the intuitive pic- 
ture which emerges from the formulation given in ((2]) of 
the Hamiltonian H: the lower is g, the more processes 
requiring the presence of no doubly occupied sites (those 
with coefficient t) are favored with respect to those re- 
quiring rid 7^ (with coefficients g and t a d)- The result 
is, instead, in contrast with what happens in the moder- 



V. CONCLUSIONS 

In summary, we have shown that the strong-coupling 
regime of Hamiltonian (pp) is described by the t-J Hamil- 
tonian in which J = 4g 2 /U. At half- filling, the effect 
of the particle-dependent tunneling rate g is to drive the 
Mott insulator from an antiferromagnet toward a param- 
agnetic configuration of fully localized atoms. At inter- 
mediate values of the interaction J7, this behavior induces 
a related behavior of the number of doubly occupied sites 
rid, which decreases down to with g, thus providing a 
criterion for inferring the effective value of g directly from 
experimental measurements of rid- 
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